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ABSTRACT 

A study on the technological feasibility of using solar concentrators for 
crystal growth and zone refining in space has been performed. Previous 
studies of space-deployed sola r concentrators were reviewed for their appli- 
cability to materials processing and a new state-of-thi-art concentrator- 
receiver radiation analysis was developed. The radiation analysis is in the 
form of a general purpose computer program. 

It was concluded from this effort that the technology for fabricating, 
orbiting and deploying large solar concentrators has been developed. It was 
also concluded that the technological feasibility of space processing materials 
in the focal region of a solar concentrator depends primarily on two factors: 

(1) the ability of a solar concentrator to provide sufficient thermal energy for 
the process and (2) the ability of a solar concentrator to provide a thermal 
environment that is conducive to the processes of interest. The analysis indi- 
cate that solar concentrators can satisfactorily provide both of these factors. 

The radiation analysis and associated computer code can determine the 
heating distribution on any target with any orientation in the focal region of a 
solar concentrator. The computer code can treat such problems as target 
placement errors, concentrator surface imperfections, etc. This computer 
code also has the capability of designing solar concentrators that give pre- 
scribed heating distributions. Faceted (Fresnel) as well as continuous sur- 
face concentrators can be treated. This analysis was used to show the 
sensitivity of the heating distribution on targets to such factors as target 
geometry, target location, concentrator surface errors, etc. It has also 
been used to illustrate that solar concentrators can be designed to give pre- 
scribed heating distributions. From these results it can be concluded that 
solar concentrators can provide the thermal environment for space processing. 
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An analysis has also been developed for determining the relation 
between mouen zone sizes and concentrator sizes. This analysis has been 
applied to three typical process materials of current interest, silicon, germ- 
anium and tungsten. 

The results of th.s study indicate that solar concentrators are attractive 
for space processing, not only from the viewpoint of providing the required 
thermal environment but also from the viewpoint of system specific power 
and system flexibility. The results also indicate that a solar concentrator 
with practical dimensions can provide the required therma 1 energy, i.e., a 
concentrator with a diameter of 10 meters (or possibly less) can process 
bars of silicon and germanium up to 10 inches in diameter. 

It is concluded that with additional analytical and experimental study, 
a solar concentrator system to provide thermal environments for all space 
processing applications can be developed. 
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Section 1 
INTRODUCTION 

Space processing offers the potential for producing some unique products 
because it provides an environment of microgravity, high vacuum and the avail- 
ability of the continuous supply of non -contaminating solar energy. The lack 
of gravity yields opportunities for redxiced, more controllable natural convec- 
tion in the liquid state, containerless processing, absence of sedimentation, 
and capillary-controlled fluid shaping and fluid motion. Reduced natural con- 
vection and containerless suspension for molten materials are highly adv?u- 
tageous in nucleation and solidification, which should lead to technologically 
advanced and economically viable materials made from metal alloys, eutectics, 
composites, solutions and pure components. The natural, high vacuum of space 
offers a contaminant-free atmosphere in which to melt, shape, mold and solidify 
process materials. The combination of the two preceding factors with the avail- 
ability of non-atmospherically attenuated solar energy provides the unique op- 
portunity for containerless melt zone processing of ultra-high melting point 
mate tls. Concentration of solar energy offers a source of high heat flux 
uncontaminated by combustion products, electric and magnetic fields, and in 
unending supply. In particular, the crystal growth and zone refining of high 
melting point materials appears very promising in space with a solar concen- 
trator configuration. 

One of the major problems of any extraterrestrial project is a sufficient 
supply of energy. For large scale operations such as the space processing 
of materials on a production basis the cost of energy in the form of batteries, 
fuel cells, etc., will be prohibitive. At this point in time the most attractive 
candidate processes for space processing involves the melting of rather large 
quantities of high melting point, high thermal conductivity, high specific heat 
materials. This will, of cou" ., require a great deal of thermal energy. It 
has been propose! to conduct the processes in the focal region of a solar 
concentrator. 
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There is economic and technical justification for pursuing the prepara- 
tion of single compound crystals (ceramic oxides and compound semiconductors) 
with maximum perfection v purity, and size to achieve high performance elec- 
tronic components and devices. Other crystals that appear attractive at this 
time are garnet crystals for magnetic bubble memories ard crystals that have 
electro-optic applications [1]. A primary advantage of zero-g for crystal 
growth is that crystals much larger than those grown on earth can be grown 
from a melt zone in space. The size of the earth-grown crystal . depend pri- 
marily on the surface tension of the material. In zero-g crystals 4 to 8 inches 
in diameter and perhaps larger may be feasible. High purity tungsten for x-ray 
targets is also an attractive product for zero-g manufacturing [2]. Tungsten 
can be zone refined using the float zone technique without crucible contamination. 

In addition to growing crystals and zone refining, some other possible 
applications of zero-g are production of biologicals, separation of materials, 
obtaining accurate values of physical properties, etc. Justification of space pro- 
cessing, candidate processes, etc., has been well studied [ 1-3] • Many of the 
processes, experiments, etc., mentioned in these studies require thermal 
energy. This of course is the function of a solar concentrator. 

The advantages and disadvantages of crystal growth and zone refining 
in space with a solar concentrator are outlined in Tables 1 and 2. Many of 
vhe points in this table have been mentioned in the text. The table is self 
explanatory. 

The objectives of this study as originally outlined were: 

• Determine the feasibility of using a system of mirrors 
to concentrate solar energy for direct processing of 
materials. 

m Define solar concentrator materials compatible with 
the space environment. 

• Develop a design for a crystal growth and zone refining 
assembly. 

• Specify materials that are most suitable to, and benefitted 
by solar processing in space. 
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Upon reviewing the literature related to solar concentrators it was found that 
the design of large solar concentrators for space applications had received a 
great deal of attention during the early 1960s. The purpose of these concen- 
trators was to power closed thermodynamic cycles for electrical power gene- 
ration. However, the technology developed for .hese concentrators is directly 
applicable to the present study. The technology develc jq included every aspect 
of the problem from materials selection and fabrication to launching and deploy- 
ment. Also upon reviewing the literature related to space processing it was 
found that a great deal of current effort is devoted to specifying materials that 
are most suitable to, and benefitted by processing in space [l-3]. Thus, some 
of the original objectives had already received much broader treatment than 
could be given to them in this limited study. Therefore the primary effort of 
this study was devoted to determining the feasibility of using a system of 
mirrors to concentrate solar energy for direct processing of materials and 
to developing a design for a crystal growth and zone refining assembly. 

To determine the feasibility of conducting zone refining and crystal 
growth in the focal region of a solar concentrator the relationship between 
the processes and the concentrator must be understood. The purpose of the 
solar concentrator if of course to provide the thermal energy for the pro- 
cesses. Upon examining the processes of interest it was found that the 
success of the processes depends to a large degree on their thermal environ- 
ment, e.g., to grow a high quality crystal, a near planar solid-liquid interface 
is mandatory. Thus the feasibility of using the solar concentrator for pro- 
cessing depends to a great extent on: (1) Can the solar concentrator supply 
sufficient thermal energy for the process?; and (2) Can the solar concentrator 
create a specific thermal environment for the process? These are the ques- 
tions primarily considered in this report. 

In the following section background information on various aspects of 
solar concentrators and the processes of interest is given. The topics covered 
are shown on the following page. 
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• Crystal growth 

• Continuous Zone Refining 

• Molten Zone Shape 

• Controlling the Thermal Environment 

• Concentrator Contamination 

• Pointing the Solar Concentrator 

• Target Placement, and 

• Previous Solar Concentrator Studies. 

Following the background material, a radiation analysis is presented 
for determining the ir radiative distribution on targets placed in the focal 
region of solar concentrators. This analysis can consider any type of axisym- 
metric concentrator and any target configuration. Also system imperfections 
such as target placement errors, optical errors, etc., are considered. 

Following the radiation analysis, a simplified heat transfer analysis is 
performed for a molten zone. This analysis allows estimates to be made of 
concentrator size versus process size. 

By using the analyses developed in this study parametric data are gene- 
rated to show the influence of various parameters on the thermal environment 
produced by solar concentrators. 
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Section 2 

BACKGROUND INFORMATION 
2. \ CRYSTAL GROWTH 

In growing crystals from a melt, crystal perfection is often a function 
ol’ the shape of the growing interface. The shape of this interface is controlled 
by the heat flow pattern in the system. The heat of fusion of the solidifying 
m aterial is liberated at the solid -liquid interface and transported away from 
the interface by conduction, convection and radiation. The relative importance 
of each of these transport mechanisms depends on the conditions, physical prop- 
erties and geometry of the system. For example, for growth above about 800 
to 1000°C, radiation from the growth interface becomes more important. This 
importance increases if the crystal has any degree of transparency in the spec- 
tral range in the vicinity of the peak of the blackbody radiation curve for the 
growth temr rature [ 4] . 

The growth interface shape is dictated by the heat flow pattern which 
is i. ont rolled by the overall thermal environment of the process. Although 
a particular heat transport mechanism is not significant at the interface, this 
mechanism may plc.y a very active role in determining the growth interface 
shape. 

The problem of the interaction of heat and mass flow with crystal struc- 
tural imperfections is quite complex and in a relatively embryonic state of 
development [4] . The relation of crystal growth to the thermal environment 
for the melt growth technique is particularly well stated by Chu [ 5] . 

•'When a crystal is grown from a melt, the solidification pro- 
cess must be controlled to avoid random nucleation. Usually, a 
single crystal seed is used, and the solidification is allowed to 
take place on the surface of the seed by adjusting the thermal 
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geometry. Under proper conditions, the growth continues the 
cry tallcgraphic orientation of the seed to form a single crystal. 
The thermal conditions at the growing liquid- solid interface are 
the most important factors governing the crystal growth process 
and are influenced by many factors. In practice, the liquid is 
heated at a temperature above the melting point by conduction 
and radiation from a thermal source. The heat of fusion is gene- 
rated at the interface at a rate directly proportional to the rate 
of solidification. Heat is lost from the interface by conduction 
through the grown crystal and by conduction and radiation from 
the crystal to thermal sinks provided by the container walls and 
any gaseous atmosphere. If there is a net accumulation of heat 
to the interface, the crystal will melt. If there is a net loss of 
heat at the interface, the crystal will grow at a rate determined 
by the rate of dissipation of heat of fusion. Thermal gradients 
in the liquid and solid determine the rate at which heat arrives 
at and leaves the interface, thus determining the rate of crystal 
growth. The control of temperature and temperature gradients 
are therefore of utmost importance in melt-growth techniques. 1 * 


Although solar concentrators have not been used to grow crystals, other 
types of thermal imaging furnaces have been used for this application [6-8]. 
Many of the findings in these studies are applicable to the present study. For 
example some of the topics mentioned are thermal environment control, oper- 
ating procedures, atmosphere control for the process, and some comments 
concerning the advantages and disadvantages of using imaging furnaces. 


2.2 CONTINUOUS ZONE REFINING 

Basically continuous zone refining is conducted by moving a molten zone 
along the length of a bar of material to bring about impurity segregation. The 
technique depends on differences in composition of the liquid and solid in equi- 
librium. By repeatedly sweeping the molten zone from one end of the bar to 
the other, the impurity concentration decreases at one end of the bar and in- 
creases at the other end. After a sufficient number of passes most of the im- 
purities are swept from one end leaving a high purity material [e.g., see Ref. 9] . 
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2.3 MOLTEN ZONE SHAPE 

An important aspect of float zone processing is the shape of the molten 
zone and the shape of the meniscus connecting the molten zone to the solid 
material. The molten zone shape and meniscus shape is dependent on the 
molten material surface tension (which is a strong function of temperature) 
and the geometry of the process, i.e., geometry of the feed stock and the pro- 
cessed material. The molten zone geometry also couples the process to the 
solar concentrator in that the geometry influences the irradiative distribution 
received, which in turn influences the solid-liquid interface, Marangoni con- 
vection. etc. Although this study does not deal with the problem of molten 
zone shape and meniscus shape the problem has received attention [ 10- 14]. 

2.4 CONTROLLING THE THERMAL ENVIRONMENT 

As mentioned previously, crystal perfection and shape depend on its 
growth thermal environment, i.e., the crystal shape is dictated by the meniscus 
shape which depends on the surface temperature, and crystal quality is a strong 
function of the solid- liquid interface shape which is dictated by the total thermal 
environment. Thus an essential part of the solar concentrator system will be 
a control system which monitors the process temperature and controls the heat 
flux distribution. Various sensing and control techniques have been considered 
by previous investigators [15-21] . 

Obtaining a prescribed irradiative distribution on a target in the focal 
region of a solar concentrator can be achieved by positioning the target and/or 
adjusting the concentrator surface. Another alternative is the use of a second- 
ary optical element to change the irradiative distribution. For example, one 
could use a ,l flux redistribute i 4 ' such as that described by Glaser, et al., [ 17] 
if uniform irradiation is desired. One could, of course, design similar second- 
ary elements to achieve other irradiative distributions. 

The amount of radiant flux available to irradiate the receiving surface 
can be controlled by use of a cylindrical shutter [ 1 6] or a Venetian blind 
arrangement [ 1 5] . 

«•. 

Convection induced by surface tension gradients. 
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Techniques are available for measuring the heat fluxes of solar furnaces 
[ 22] and the temperature of targets in the focal region [ 21 , 23] . These meas- 
uring techniques can be coupled to controls for achieving the desired thermal 
environment. 

2.5 CONCENTRATOR CONTAMINATION 

There are many materials available for concentrator reflective surfaces 
[ 24-27] . However it seems that none would be immune to degradation due to 
protons and ultraviolet radiation [ 27, 28] . A protective coating would be essen- 
tial for long term surface optical quality. The reflective sur^ce can also be 
degraded by contamination due to process material outgassing. This can be 
prevented by shielding the reflector by enclosing the process area with glass, 
or if this contamination is very slow, it could be removed periodically using 
RF -excited oxygen plasma [ 10] . Analytical models are available for com- 
puting the adsorption and desorption of materials under specified conditions 
[29,31]. 

It has alto been suggested that a gas jet could be used to deflect out- 
gassed molecules to reduce contamination if the process is to be conducted 
in the vacuum of space [32]. 

2.6 POINTING THE SOLAR CONCENTRATOR 

It has been stated that a pointing accuracy of about 1 minute of arc is 
needed for a high quality solar concentrator [ 32] . Several devices have been 
developed for pointing solar furnaces [ 33-35] and appendages to space vehicles 
in general [ 36] . The pointing accuracy requirements for a solar concentrator 
arc much more lenient than those for telescopes. Thus it seems that there 
will be no problem in achieving the required accuracy. However, the accuracy 
is such that it may be necessary to dynamically isolate the solar concentrator 
system from any manned laboratory. 
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2.7 TARGET PLACEMENT 

The spatial location of the process or target in the focal region of the 
solar concentrator is of extreme importance. The rays being reflected onto 
the target are converging on the target from all portions of a solid angle that 
will probably exceed 4 n steradians. The n focal region 11 of the concentrator 
is a relatively small volume. Thus any rotation or translation of the target 
from the desired position can cause great variations in the heating distribu- 
tion and the magnitude of the heating. 

2.8 PREVIOUS SOLAR CONCENTRATOR DEVELOPMENT STUDIES 

Several efforts have been directed toward the development of solar con- 
centrators for space applications f 24, 37-44] . These efforts have been con- 
cerned with many aspects of the problem. They have covered fabrication 
techniques, reflective coatings, reflective surface degradation, deployment, 
surface accuracy requirements, etc. 

Methods of fabricating and testing concentrators of practically any size 
seem to be well developed. A variety of deployable designs have been con- 
ceived, e.g., petal configurations, Fresnel type reflectors, umbrella configu- 
rations, and inflatable types [ 39,41, 45] . Most of the compact designs for 
deployment have optical accuracies sufficient for obtaining concentrations that 
that would be suitable for many applications. However, the one piece non- 
folding type concentrator has the advantage tha 4 ; ts performance can be well 
characterized before it is launched. 

Electroformed nickel concentrators appear to the most efficient one- 
piece concentrators. However, stretch-formed aluminum has a higher 
specific power based on weight. For a given power requirement an aluminum 
concentrator would have a greater diameter [41]. 
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Rather comprehensive studies of solar concentrators for orbiting were 
conducted in Refs. [43] and [44] . In Ref. [43] the Air Forced EROS Program 
(Experimental Reflector Orbital Shot) is described in which a 4.5 foot diameter 
reflector of the Fresnel type was successfully orbited. The effort by Stewart 
[44] was concerned with paraboloid concentrators up to 30 feet in diameter. 
This is a rather comprehensive study that covers many aspects of solar con- 
centrators that are of interest in this study. Although a complete review of 
this work cannot be given here some points of particular interest will be listed. 


• The concentrator has a specific weight of 0.87 pounds per 
square foot of projected area* 

• The reflective coating is chemically deposited silver with 
a reflectivity of at least 0.91. 

• Blade grinding or plastic overlay master fabrication tech- 
niques are best for obtaining highly accurate, highly reflec- 
tive concentrators. 

• Calculated space environment degradation of the reflective 
surface will be less than 1% in the reflectance in 10,000 
hours. 

• The concentrator developed was ground tested and has good 
launch and orbit capability. 
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Section 3 

RADIATION ANALYSIS 

Analyses for determining the irradiation on receivers placed in the 
focal region of solar concentrators have been performed by many investiga- 
tors [46-57J . If the attributes of all of these analyses were combined one 
would no d'. ubt have an analysis of sufficient generality to treat any solar 
concentrator -receiver combination with consideration given to the various 
details such as limb darkening effects, optical errors, positioning errors, 
etc. Some of these analyses treat only flat receiving surfaces located in the 
focal plane of a paraboloid. Some do not consider variation in the solar disc 
radiance (limb darkening), most consider only paraboloid concentrators, etc. 
Thus, even though all aspects of the problem have been treated, it seemed 
that a general unified analysis did not exist. Thus an analysis was developed 
that is very general and considers as many aspects of the problems as possible 
and is applicable to any solar concentrator problem including the geometries 
that are of interest to crystal growth and zone refining. A computer code has 
been developed for this analysis. 

The analysis consist primarily of developing ray trace expressions and 
energy balance relations that take into consideration the characteristics of the 
source. The radiation analysis is the heart of the computer code for solar 
concentrator analysis and design. A general treatment will be made of the 
radiation analysis. The analysis will then be modified as needed to treat the 
problem at hand. 

The radiation from the solar source that is important to the thermal 
image in a solar concentrator is in a relatively narrow bandwidth (0.3/ji to 
3 fJL). Thus in the analysis that follows, the spectral properties of the radia- 
tion will not be considered. 
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3.1 RAY TRACE 

Consider the reflecting surface described by the implicit function S , 
and the receiving surface S£ as shown in Fig.l . The law of reflection re- 
quires that Pj = P| and that Vj , and Nj are coplanar. These conditions 
can be stated mathematically as follows: 





N, = 0 


( 1 ) 


( 2 ) 


If the ray vector V ^ originates at a point 'O' then position vectors for the ray 
vector extremities can be written as follows: 


R =ix +jy +kz 
o o J 'o o 

(3) 

Rj =^ x ! + j Yi + k Zj 

(4) 

^2 =Tx 2 + "^ y 2 + irz 2 

(5) 

With the above the ray vectors can be expressed as: 


V = R, - R 
1 1 o 

(6) 

W 

(7) 
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This of course gives the ray vectors in terms of the coordinates of the source, 
reflecting surface, and receiving surface. B"- making use of the fact that the 
gradient of a scalar function describing a surface is a vector normal to the 
surface, the ray trace equations can be expressed in terms of ray coordinates 
and partial derivatives of the reflecting surface, i.e., Nj or its negative image 
(depending on the function S ^ ) is given by 

Nj = V S, (8) 

where 




(9) 


With the relations above it can be seen that the ray trace equations (1) and (2) 
are partial differential equations in terms of the source, reflector, and receiver 
coordinates and the reflector surface derivatives. Equations (1) and (2) can be 
expressed in functional form as follows: 


g(R o , Rj, R 2 , VSj) - 0 


HO) 


h(R , R. , R, TS, ) = 0 
o 1 ? 1 


(ID 


The above describes the path of a ray through any reflecting optical system. 
The application of these equations to the problem of interest will be made 
later. 


3.2 ENERGY BALANCE AND IMAGING THE SOLAR SOURCE 

Consider the arrangement in Fig. 1 where S q is the radiation source, 
Sj is the reflecting element and S 2 is the receiving surface. There, are two 
ways to view the imaging of the source S q onto the receiver S 2 through the 
reflecting element S^. Let us explore these. 
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First, consider an emitting patch on the source that is sufficiently small 
relative to the characteristic dimensions of the system to be considered as a 
point source (see Fig. 2). The image of the reflector that is produced on the 
receiving surface depends on the geometry of the reflector (including curva- 
ture, etc.) and the geometry of the receiving surface on which the image is 
formed, and of course the relative spatial location of the source, reflector, 
and receiver 

If this patch has area uA q and radiant emittance M then the radiant flux 
from this area will be 


4> - MdA 

o o 


( 12 ) 


Consider an element of the reflecting surface dA j located relative to the source 
as given by N q , Vj, and . The radiant intensity of the source can be expressed 
as a function of 0 , x.e. 


I(f> 0 >=I n *(f> 0 ) (13) 

where 1 is the intensity in the direction of N . For a prescribed diB ), I 
n o o n 

can be determined from the total radiant flux by integrating over the half space 
seen by dA Q , i.e. 


or 



(14) 


1 

n 


7 


4 > 




(15) 
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Thus the radiant flux intercepted by dAj will be 


*1 * 


where the solid angle dft| is given by 

dA, cosP, 

jo * *■ 



or with (12), (15) and (17), (16) can be written 


( 16 ) 


( 17 ) 


MdA + (P ) cosP.dA. 

<». = 2 - . i L (18) 

/#0 0 )dfl Vj-Vj 

4r 

Thus if the source is characterised, (18) can be used to compute the radiant 
flux on dA j due to the emitting element dA Q . 

The radiant flux available to irradiate dA^ is after being attenuated 

by dAj . If dA j has reflectivity >'j, which can be dependent on Pj, then 


*2 = y l (P 1 } *1 


(19) 


The irradiation on the receiving surface can be determined by differentiating 
the ray-trace equations (10) and (11) and coupling them to (19). If the source 
is extended then the irradiation at a point on the receiver is determined by 
integrating over the source through the reflecting surface. One may *nt v 
this as the overlaying of reflector images produced by point sources as shown 
in Fig. 3. 


18 


LOCKHEED HUNTSVILLE RESEARCH « ENGINEERING CENTER 



Reflector 


LMSC-HREC TR D390666 



19 


LOCKHEED ■ HUNTSVILLE RESEARCH & ENGINEERING CENTER 


Fig. 3 - Overlaying Images Produced by Point Sources 



LMSC-HREC TR D3?0666 


Let us now consider an extended source being imaged onto a receiving 
surface. We will assume that the reflecting element is sufficiently small 
relative to the system characterizing dimensions that it will act as a 'pinhole,' 
i.e., the image formed on the receiving surface will b< determined by the 
source shape aad the receiving surface rather than by the shape of the reflec- 
ting surface. This situation is depicted in Fig. 4. Hie image geometry for 
each 'pinhole' is determined from ray trace equations and the geometry of 
the receiving surface. However, the irradiation within the image is dependent 
on the source characteristics as well as geometric considerations. Thus the 
irradiation at a point on the receiving surface is determined by integrating over 
the reflector while considering the source characteristics. One may view this 
as the overlaying of images produced by the source geometry as shown in Fig. 5. 

The foregoing can be summarised as follows: the radiant flux incident on 
dAj can be computed from (18) by integrating over the source S^, and that avail- 
able to irradiate the receiver can be obtained by integrating (19). For the case 
of a point source the ray trace equations must be differentiated holding the 
source coordinate constant. This gives the mapping of the reflector onto the 
receiver. By coupling this to energy considerations an expression for the ir- 
radiation can be derived. However, for the situation of a pinhole reflector the 
reflector coordinates are held and the ray trace equations are differentiated. 
This gives the mapping of the source onto the receiver. By coupling this to the 
characteristics of the source an expression for the irradiation can be derived. 

Thus there are basically two approaches that can be used to compute the 
thermal image of a radiation source such as the sun. Due to the characteristics 
of the source and its size relative to its distance from the reflecting surface, 
the second approach is more desirable for this application. 

3.3 THE SOLAR SOURCE 

Considerable effort has been devoted to characterizing the solar source 
in terms of its irradiance outside earth's atmosphere (solar constant), spectral 
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properties, and variation in radiance across the source [50, 58-6i]. The 
most generally accepted value of the solar constant is 1.94 calories per square 
centimeter per minute. The thermal radiation from the solar source is pri- 
marily in the wavelength range of 0.3 to 3.0 microns. This is a relatively 
small wavelength variation. As a result the spectral properties ot the radia- 
tion can normally be averaged for the purpose of computing a thermal image. 
However, the brightness or radiauce as a function of position on the solar disc 
is very important in determining the irradiation distribution on targets in the 
focal region of solar concentrators [50]. There are several models of the 
solar disc available [50, 58]. The solar disc is normally assumed to subtend 
an angle of 32 minutes when viewed in the neighborhood! of earth. However, 
there are slight seasonal variations in the solar disc size due to the change 
in earth-sun distance. This variation is from 32.59 minutes in January to 
31.52 minutes in July [57]. The effective temperature of the solar source 
as seen from our planet is somewhat controversial, however an often quoted 
number is 5800°K. 

To describe the solar source let us consider a very small reflecting 
element imaging the solar disc onto a plane normal to the reflected ray vector 
as shown in Fig. 6. The radiance of the sun decreases from its center to the 
outer edee and thus the irradiation in the solar disc image varies in the same 
manner. Expressions for the intensity within the sun image integrated over 
all wavelengths have been given by Kamada [ 50] and Angstorm and Angstrom 
[58]. 


The variation in intensity can be expressed as a function of the angle 
between the ray vector of interest and the image centerline ray vector. 
Letting any ray vector to the sun image be the intensity distribution 
across the sun image will dec rease radially as shown in Fig. 7. The in- 
tensity within f he sun image is expressed in terms of the image center 
intensity and the angle Intensity is defined in Eq. (20) shown on the 
following page. 
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Fig. 6 - Imaging the Solar Source 
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! 


i 



(20) 


For an element of the image dA^ located at the termination point of V c 
(assuming £ is small) 


dfi = dA 2 /(V^.V^) 


( 21 ) 


Irradiation is defined by 


E 



(22) 


Thus with the above *he irradiation on the surface is related to the intensity 
by 


or 


E(|) = I(|)/(V^.V^) 


(23) 


E(|) = E o G(fc) 


The irradiation on the plane being considered must satisfy the energy balance 


-J 


EdA 


(24) 


The distance from the center ray to any point on the image can be expressed 


v 



(25) 
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Now 


dA = Zn udv 


2 7r£ d£ 



(26) 


Thus 


F rom which 


d *2 


2ttE 

o 


V -V 
2 2 



y* G(g)£d| 


E 

o 


7 

2 arV 2 *V 2 j G(§)£d£ 


(27) 


(28) 


Thus with a prescribed intensity distribution and a prescribed reflecting ele- 
ment which gives the reflected radiant flux d$, one can compute the irradia- 
tion at any point in the solar image produced by the reflecting element with 
the use of (23) and (28), i.e.. 


dt ? G(|) 

E(|) = r (29) 

_ _ f° 

V 2* V 2 J G <$>$ d & 

o 

The above gives the irradiation distribution on a plane normal to the 
reflected ray vector. This will be extended to non-planar receivers in another 
section. Equation (29) characterizes the solar source imaged by a small 
reflecting element. 
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The effects of limb darkening are important u nigh concentrations are 
achieved. In situations where the radiation is spread out over several sun 
image diameters, the limb darkening effects are obliterated, 

3.4 APPLICATION OF THE ANALYSIS TO SOLAR CONCENTRATORS 

The analysis developed in the previous sections will be applied to deter- 
mine the heat flux on targets in the focal region of solar concentrators. The 
analysis developed is rather general and is applicable to any solar concentrator 
situation. However, for the particular situation under consideration the general 
analysis is not needed. An abbreviated form can be used without any loss of 
genoraL;,. The concentrator and receiver configurations of interest are axi- 
syrimetric. However, non-axisymmetric conditions will exist in the system 
due to pointing errors, optical imperfections, and target placement errors. 
However, for the particular situation at hand, a general analysis is best accom- 
plished by using an axisymmetric analysis that is modified to consider the non- 
axisymmetric conditions. 

For the axisymmetric configuration in Fig. 8, the coplanarity condition 
given by Eq. (2) will automatically be satisfied since all rays are meridional. 

Let the coordinate system be oriented such that the unit vector k is pointing 
toward the solar source. This of course implies 

V L = - k (30) 

if s and r are taken as the concentrator and receiver axisymmetric coordi- 
nates, respectively, and e as a unit vector in the plane of symmetry (see Fig. 8) 
then the position vectors of the ray extremities given by (4) and (5) become 


Rj = e s t k Zj 


( 31 ) 
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- e r + k 


(32) 


and the reflected ray vector ( * ) becomes 


V 2 = e (r - s) + k(z 2 - Zj) 


(33) 


For the axisymmetric situation 


s = S (Zj) 


(34) 


cr the implicit surface function is 


Sj (s, Zj) = s - s (Zj) = 0 


(35) 


and 


Nj = TSj = e - k s* 


(36) 


Thus with (36), (33), a 4 \d (30). the ray trace equation (1) becomes 


-k 


e (r - s) + k (z 2 - Zj j 

>/(r - s) 2 + (z 2 - Zj ) 2 J 


(e - k s') = 0 


(37, 


or upon performing the dot product 


ll 


s’ || V 2 -(z 2 -Zj) | r ir - s) = 0 


( 38 ) 
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where 




s r + (2, 


V 


(39) 


Equation (38) gives the ray path for source centerline rays. Other rays 
from the source are described by the above and the source geometry. 

Let us relate the axisymxnetric coordinate s to the coordinates y and 
z. Defining 


4# = 

-1 z — 

cos ( e • j ) 

(40) 

y i = 

S COS' ' 

(41) 

x i = 

s si riy 

(42) 


Con« :J .er a differential element of the reflecting irface as shown in Fig. 9 
the area of this element can be expressed 


dA 


1 


s [ (ds) 2 + (dzj) 2 ] 



(43) 


or 


dA = s[(s*r + lj 


1/2 


dy dz 


1 


If the irradiation on „ plane normal to V ^ is given by the solar constant K 
then the radiant flux intercepted by dA^ will be 

d? j = KdAj cosjij (44) 
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where Pj can be defined by 



By combining (43). (44), and (45) 

d$j = Kss* d+ dSj 


(45) 


(46) 


If the reflecting element has reflectivity which can be a function of , 
the radiant fhuc reflected onto the receiver will be 


d4> ^ = <|)j (Pj) Kss'diji dZj 


(47) 


It was shown previously how the irradiation on a receiving surface 
normal to the reflected ray could be computed if the incident radiant flux 
d $2 and the source, description G(|) were known. Let us now consider a 
general receiving surface onto which the solar source is being imaged as 
shown in Fig. 1C. We wish to compute the irradiation at an arbitrary point 
p. Looking at a plane containing 0 and p as shown in Fig. 11, the irradiation 
on a plane normal to and passing through p can be computed. If the co- 
ordinates of p are (x , y , z ), V, can be written as 
K P P P I 



(48) 


where 



ix +jy + k z 

P J r P P 


i s sirn^ + j s cos4* 


+ k z 


(49) 

(50) 
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With the above £ can be computed from 


£ = sin* * 


Fz - %\ 


(51) 


The angle P q between the incident ray and the surface being irradiated 
can be computed from 


P 



(52) 


where 


N 2 = vS 2 (x p'VV 


The irradiation at p is given by 


Ep = E(£) cosp 2 


(53) 


(54) 


With (29) and (47) the above becomes 


y j (P j ) K s s' G(£ ) (cos P 2^^ i 

E - T 

p i. 


Zn V t -V £ 


/° 

/ G(£ ) £ d£ 


(55) 


The above gives the irradiation on an arbitrary receiving point with any orien- 
tation. This of course is for a differential element of the reflecting surface. 
To determine the total irradiation at a point the above must be integrated over 
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and 7.j 


or in other words over the concentrator surface S 


1* 


Let 


Y. (P ,) Kss'G(|)co8p 2 

T = - 1 1 (56) 

Z* V| . W % f 

o 


Then 


"tot 


=■£ % = //ms. 0 2 > Td+d«! 


y z. 


(57) 


where 6(| . P^) is defined as follows 


M£.p 2 > = 


1^ 

I i */2 


0 


(^ > ^o 

| P 2 > r/2 


(58) 


£ and are dependent on ^ and 7j. The purpose of 6 (| , P^) is to take care 
of two points that must be considered in the integration process. These are: 
(1) if the receiver point lies outside the image formed by the reflecting ele- 
ment there is no contribution to the irradiation and (2) if the receiver point 
cannot 1 see 1 the reflecting element there is no contribution to the irradiation. 

3.5 TARGET PLACEMENT ERRORS 

The ir radiative distribution on a receiving surface placed in the focal 
region of a solar concentrator can be changed significantly by small changes 
in the spatial location and orientation of the receiving surface. If the receiver 
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is initially in the axi symmetric location any placement errors or intentional 
movements of the receiver can be described by a rotation in one plane and 
a translation along the axis normal to this plane as shown in Fig. 12. 1 he 

new coordinates and orientation of the receiver can Le expressed in terms 
of the original coordinates and surface normal vectors as follows. 

The surface normal vector is expressed 

^ ^ dS ? 9S ? aS ? 

N 2 = ^~5F + j Tf + * "5T (59) 


where S£(x, y, z) is the implicit function describing the receiving surface. 
The component of N, in the x direction will not be changed. Thus let us con- 
sider only the components in the v-z plane being rotated as shown in Fig. 13. 
We can express the components in the y-z plane as 


n = j 


as 2 

W 


as. 


+ k 


dz 


(60) 


Let us express n after rotation as 



(61) 


Nowj n'^l n| and we can express the angles in Fig. 13 as follows 


DS-, 

cosfi = 



( 62 ) 
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and 



By making use of the identities 


cos(M + a) = cos/JL cosO - sin fl sinO) 
sin ( /LX + u.) = sin |i cosu) + cosflsino) 


and the above we can show that 


/asA* as. as 7 

(~3 — ) = -T COSUu - — 5 — si 

\ oz f dz a y 




COSU) 


aS 2 

~9z 


SUIju 


Thus the new surface normal vector is 



cosa) + 


as. 


ciZ 


sinai^ 


+ k 


as 2 

9z 


CCSJJ 


3y 


sinto\ 


( 66 ) 


(67) 


( 68 ) 
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and the new coordinates of a point on the surface can be expressed 


x* = x + A 

y* = y cosu) + z sinu) 

z* - z cosoo - y sino) (69) 


3.6 OPTICAL ERRORS 

Any reflecting surface will have some degree of geometric imperfection 
or waviness that occurs because of fabrication techniques, stresses, etc. The 
exact nature of these imperfections are very difficult to determine a priori. 
These imperfections do not normally follow a normal statistical distribution 
and are probably very dependent on the individual reflector. It has been 
pointed out by Stewart [44] that the imperfections cannot be simulated by a 
Gaussian distribution. However, if data describing the geometric imperfec- 
tions is not available the use of random surface errors as described by a 
Gaussian distribution can give some insight to the influence of geometric im- 
perfections. The approach for randomly introducing surface errors is illus- 
trated in Fig. 14. The surface slope is perturbed randomly between 0 and - 
three standard deviations in a Gaussian manner. The angular imperfections 
corresponding to one standard deviation must of coarse be specified. 

3.7 REFLECTOR DE r .GN 

By coupling the geometric optics and energy balance analyses to finite 
difference formulas, reflectors that give prescribed ir radiative distributions 
can be designed. The idea of overlaying the thermal images produced by 
differential elements of reflecting surface is utilized. If an irradiative distri- 
bution is prescribed, a reflector surface is determined by *building up 1 the 
reflector by forward differencing schemes. The irradiative distribution 
requirements and source characteristics guide the differe* ^ing scheme in 
arriving at the surface shape. This idea is illustrated in Fig. 15. 
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While this technique can be applied to any source-receiver combination 
it should be pointed out that there are inherent limitations when an extended 
source is used. If a point or collimated source is used the reflector is de- 
fined by a differential equation as shown in [62-63] and the only constraining 
factor is the energy available. However for an extended source an integral 
equation defines the reflector shape. For this situation the source character- 
istics determine whether a particular ir radiative distribution can be achieved. 
For exctinple, as pointed out in Ref. [63] irradiative distributions with la*ge 
gradients along the receiver surface may be very difficult or impossible to 
achieve. 

3.8 PARAMETER DEFINITIONS AND NON -DIMENSIONA LIZING THE 

EQUATIONS 

It is rather easy to non-dimensionalize the equations developed in the 
analysis for solar concentrators. Only a characteristic dimension and a 
characteristic radiant quantity are needed. The obvious choice for the char- 
acteristic radiant quantity is the solar constant K. The characteristics dimen- 
sion chosen will be a focal length f. Although this quantity is somewhat neb- 
ulous for concentrators without a true r focal point* it is the quantity that will 
be used. Most of the results published in the open literature are for para- 
boloids. An important parameter for any optical device is the aperture ratio 
— which is the ratio of the aperture or opening that receives radiation to the 
focal length. Results are normally presented as a function of this parameter. 

Hiester, et al. t [47] defined the concentration ratio as the ratio of the 
heat flux at the sun imagj of a sun furnace to the actual flux received from 
the sun ac normal incidence after reflection. This definition can be applied 
to any point on a receiving surface. Actually this definition will become 
evident up~m non-dimensionalizing the irradiation given by Eq. (57) i.e. f 



( 70 ) 
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Another useful parameter defined in [47] is the concentration efficiency. 

This is the ratio cf the power reflected by the concentrator to that received 
by the target. That reflected can be obtained by integrating Eq. (47), i.e.. 


* „ - 


//*> // 


yjiPjjKss* «*)■ d*j 




»»» z. 


(71) 


The power incident on the receiver is obtained by integrating the irradiative 
distribution on the receiver, i.e.. 


$ 


r 


•// 

z 6 


E (z,e ) dzd6 


and the concentration efficiency is given by 


(72) 


1 = 


(73) 


This of course takes into consideration radiation 'spillage* due to target 
geometry and orientation. 

A parameter that describes the amount of radiation actually absorbed 
will not be defined. Such a parameter would depend on too many factors 
concerning the target surface properties and shape and would not be good 
for making comparisons in general. 

A catch-all factor called the 'furnace factor' is often defined for a situa- 
tions such as that being considered. It is felt that the definitions of this quantity 
in the literature are not consistent, therefore it will be avoided. 
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Section 4 

PROCESS HEAT BALANCE 

The solar concentrator analysis computer code simulates the concentrator- 
receiver configuration and predicts the heat flux distribution on the molten zone. 
However, the geometry of the molten zone, the optical properties of the process 
material surfaces, and the thermal and physical properties of the process ma- 
terial dictates the thermal energy requirements and thus the solar concentrator 
size. Therefore a heat balance on the molten zone and the connecting material 
has been performed. 

It has been mentioned previously that c rystal quality and shape depend 
on meniscus shape and the solid-liquid interface shape. The meniscus shape 
depends on geometry and surface temperature. The solid-liquid interface 
shape, however, is very implicit in the total thermal environment of the pro- 
cess. The detailed thermal analysis to predict the solid-liquid interface is 
rather formidable. Such an analyses must be coupled with the solar concen- 
trator analysis to optimize solar concentrator shape, positioning of the pro- 
cess, etc. However, such an analysis did not fit within the time frame of the 
present study. Therefore the simplified one -dimensional analysis presented 
in Ref. [9] will be modified for estimating concentrator size versus process 
size. 


If the isotherms in a cylindrical ingot are very nearly flat and perpen- 
dicular to the axis of symmetry, the following one -dimensional differential 
equation results for the heat balance (see Fig. 16) 


K 



- VpC 

P 


dT 
s 

dz 



( 74 ) 
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where 

T s = surface temperature 
7 . = axial coordinate 

k = conductivity 

V = velocity of the solid -liquid interface 
p = density of the material 
Cp - specific heat of the material 
P = perimeter of the bar 
A = cross-sectional area of the bar 
q = net heat flux at the surface 


If the solar concentrator is directing radiation only to the molten zone 
and if the process is being conducted in the vacuum of space, the heat losses 
along the solid bar are primarily due to radiation. This can be expressed as 


q = oc F(T s 4 - T a 4 ) (75) 

whe re 

o = Stefan -Boltzmann constant 

€ = emissivity of the material 

F = radiation geometry factor 

T = ambient temperature 

cL 


For the situation under consideration, T << T , thus 

a s 

q = ot FT 4 
^ s 


l?6) 


If the melt 
P/A = 2/R 


velocity is very small, with the aid >f Eq. (76) and noting that 
for a cylinder, Eq.(74) becomes 



(77) 
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The solution to this differential equation is 



whe re 


Tj = freezing temperature of the material 
R - radius of the cylinder 

Thus the heat loss along a cylinder of length L can be expressed 


(78) 


Q =2*Rcof/t (z)dz 
s J s 


(79) 


and by using Eq. (78) and integrating 


Q = Re oFT? W 

S D f 


1 - 


(1 + L/W) 


573 


(80) 


where 



(81) 


If appreciable stirring occurs in the molten zone (this will probably be 
the case since Marangoni convection will be present) then the surface of the 
molten zone can be assumed to be at the freezing temperature . The radia- 
tion from the molten zone surface will be 


;/ 


Q m = c oFT f / p(z)dz 


( 82 ) 
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where 


L = molten zone length 

P(z) = perimeter of the molten zone as a function of z 


The integral is the surface area of the molten zone, thus let us define 


A m = f P(z)dz 


(83) 


Thus with the aid of Eqs. (81), (82) and (83), the total heat loss by the process 
will be 


Q t = Q m +Q , + Q ? = Of F T/A^ 

t m si m m f m 


♦?» R i £ i° F i T f w i 


i - 


(1 + L/Wj) 


573 


*!* R fz° F 2 T t W 2 


(I 4 V W 2» 


m 


(84) 


where 


W l = 


5R l k l 


1/2 


9C l° F l T f 


(85) 


W 2 


5R 2 k 2 
9€ l° F 2 T f 


1/2 


( 86 ) 


With the above analysis the thermal energy requirements for a particular 
molten zone can be established and the solar concentrator can be sized. 
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This analysis offers some flexibility in that a transition in bar size can be 
treated and the molten zone can be of arbitrary shape. 
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Section 5 
RESULTS 

With the analytical tools developed in the previous sections several 
sets of computations have been made to illustrate the influence of various 
factors on the thermal image produced on a receiving surface by a solar 
concentrator. The factors considered are: 

a Target Geometry 

• Target Location 

• Concentrator Size 

• Concentrator Types, and 

• Optical Errors 

Before computations can be made a model of the solar source must be 
specified. As mentioned previously, several models have been proposed; 
however, for the calculations made here the model proposed by Kannada [ 50] 
will be used. This model is described in Fig. 17. The solar constant for 



P ig. 17 - Solar Disc Model 
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2 

these calculations will be 1.94 cal/cm -min. The results are presented in 
terms of concentration ratio. This allows the results to be easily applied 
for situations where the solar constant is not 1.94 f e.g., terrestrial applica- 
tions. The concentration ratio can be converted to watts per square centi- 
meter by multiplying by 0.1353. 

When radiation is used for heating, computation of temperatures that 
can be achieved in general is meaningless. The characteristics of materials 
that determine the temperatures attainable, temperature gradients that will- 
be established, etc., are both thermal Mid optical [ 1 8] . For a prescribed 
situation, the net heat flux into a mate, al depends on its surface properties, 
i.e., absorptivity and reflectivity. The internal temperature distributions 
depend on the specific heat and thermal conductivity. Thus the results pre- 
sented will not be in terms of temperature. Irradiation distributions will be 
presented in terms of concentration ratio. 

The simplified one -dimensional heat transfer analysis in Section 4 has 
been used to approximate bar sizes that can be melted for prescribed situa- 
tions. Although this analysis is not exact, the results presented are good 
approximations. These calculations have been carried out for tungsten, 
silicon and germanium. 


A comparison of Fresnel concentrators with continuous surface concen- 
trators is made. Also a concentrator is designed which gives a near uniform 
irradiation distribution. 

The results presented are in non-dimensional form. As mentioned pre- 
viously, all spatial variables and parameters have been non-dimensionalized 
by the focal length f. Thus z= 1.0 corresponds to the focal point. The re- 
ceiver dimensions have also been divided by f. Thus the results are for any 
geometrically similar system. 
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S.l THERMAL IMAGE ON TARGETS NORMAL AND PARALLEL TO THE 
OPTICAL AXIS 

To develop some insight into the radiation environment in the focal 
region of a solar concentrator the concentration distribution on targets 
normal and parallel to the optical axis has been computed. These targets 
are located at different distances from the optical axis and from the focal 
point. T^ make comparisons with previous work, e.g.. the results of Kamada 
[50] it was necessary to consider a paraboloid concentrator. Thus the curves 
in Figs. 18 and 19 are for a paraboloid with a specular surfac and an aperture 
ratio of 4 (i.e.. a = 4). 

In Fig. 18 the concentration distribution on cylinders of various radii is 
shown. The magnitude of the concentration decreases with increasing radius. 
This, of course, is due to the surface being irradiated being moved further 
from the focal point. If the cylinder centerline is located on the optical axis 
and has a radius sufficient to intercept all of the radiation, then one may view 
this decrease in concentration with increasing radius as a result of spreading 
the same amount of energy over larger and larger areas. 

In Fig. 19 the concentration distribution is given for a flat target placed 
in the focal region of the paraboloid normal to the optical axis. For this situ- 
ation the concentration distribution is symmetrical about the optical axis. The 
parameter u in the figure gives the location of the receiving surface relative 
to the focal point. In non-dimensional form 

u = z - 1 . 0 

Thus, for values of u< 0 the receiver is between the focal point and the con- 
centrator. The concentration distribution of course has the greatest value 
for u = 0. The curves for u/ 0 may be thought of as !, out of focus" conditions. 
Note that the distribution magnitude is somewhat smaller for positive values 
of u than it is for the negative values with the same absolute value. This 
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Fig. 18 - Concentration on a Right Circular Cylinder Located on the Optical 
Axis of a Paraboloid with an Aperture Ratio of 4 
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Fig. 19 - Concentration on a Flat Receiver Normal to the Optical Axis 
of a Paraboloid of Aperture Ratio 4 
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asymmetry is due tc the greater distance over which the sun rays have to 
diverge before striking the receiver. 

Although the results in Figs. 18 and 19 cannot be used to determine the 
concentration on any surface placed in the focal region, they qo give some 
insight to the radiation environment in the focal region. One can see that 
the positioning of a surface relative to the focal point and surface orientation 
can drastically influence the irradiative distribution on the surface. Thus 
the process geometry and heating distribution are strongly coupled. 

5.2 THE INFLUENCE OF PLACEMENT ERRORS ON THE THERMAL 
IMAGE 

In the previous section the influence of target location on the thermal 
image was illustrated for axisymmetrical situations. The influence of asym- 
metrical target placement will now be considered. Again the paraboloid con- 
centrator with an aperture ratio of tour will be used. To make meaningful 
comparisons t e target geometry must be held constant. Since a cylindrical 
receiver is of interest for space processing applications, a right circular 
cylinder with R = 0.01 will be used. This radius was chosen because it is 
near the lower limit of cylinder sizes that can intercept all the radiation 
(i.e., rj = 100%). (See Fig. 20). 

Two types of placement errors can occur as indicated in the analysis 
in Section 3.6. Rotation of the target will be considered first. The center 
of rotation will be assumed as being the vertex of the paraboloid, i.e., 
z = 0. Small rotations from 2 to 8 minutes have been considered. These 
were chosen because it was felt that they are realistic for the systems being 
considered. Note in Figs. 21 t 1 ^ough 24 that the co. entration distribution 
amplitudes begin to vary mc/e . mnd the receiver as j increases. The 
concentration at the 0=0 position decreases. This side of the cylinder is 
moving away from the focal point. The 0= n side of the receiver is moved 
closer to the focal point as the cylinder is rotated. With rotation only the 
heating distribution for z < 6 < Zr is identical to that for 0 <_ 6 < tt. 
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Fig. 20 - Concentration Efficiency for a Paraboloid with Aperture Ratio 
= 4 Irradiating a Cylinder 
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Fig. 22 - Concentration on a Cylindrical Target Rotated 4 Minutes 
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The concentration distributions in Figs. 25 through 28 are the result 
of translating the receiver from the optical axis from 0.002 to 0.008 units. 

The results are similar to those obtained when the receiver was rotated. 

The variation in the concentration amplitude is again due to the position of 
the various sides of the receiver relative to the focal point. The translation 
is along the x-axis. Thus, the heating distributions are identical f}r 3x/2< 0 
< 3r/2 and ^/2< 0<_3?r/2 

Translation and rotation of the target can be treated together. For 
this situation there is no symmetry in the heating around the cylinder. To 
illustrate, results for this type situation are given in Fig. 29- 

Thus the heating distribution on the process material is strongly in- 
fluenced by any inaccuracies in process placement in addition to process 
geometry. 

5.3 OPTICAL ERRORS 

The influence of imperfections in the geometry of the concentrator has 
also been considered. As pointed out earlier, the surface errors do not 
normally follow a normal distribution. However, to illustrate the influence 
of concentrator surface errors. A normal distribution of random imperfec- 
tions has been used. Figure 30 illustrates the influence of various amounts 
of surface error on concentration distribution; o is one standard deviation 
of the surface in minutes of arc. As one would suspect these surface errors 
tend to ’’flatten" the distribution. 

5.4 A COMPARISON OF FRESNEL AND PARABOLOID CONCENTRATORS 

Fresnel concentrators have the disadvantages that all o: their area 
cannot be utilized and they cannot produce concentrations as high as an equiv- 
alent paraboloid. The loss in usable area is inherent in the design as shown 
in Fig. 31. With the same amount of usable reflecting area a Fresnel concen- 
trator cannot produce concentrations as high as a paraboloid. This is due to 
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Fig. 25 - Concentration on a Cylindrical Target T ranslated 0.002 I'nits 
from the Optical Axis 
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Fig. 26 - Concentration on a Cylindrical Target Translated 0.004 Units 
from the Optical Axis 
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Concentration on a Cylindrical Target Translated 0.008 Units 
from the Optical Axis 
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Fig. 30 - Influence of Optical Errors on Concentration Distribution 
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Fig. 31 - Area ’utilization for Fresnel Concentrators 
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the reflecting elements being farther from the focal point. The sun rays have 
a greater distance over which to diverge. However the Fresnel geometry is 
desirable in some situations* e.g.. a foldable concentrator. The results in 
Fig. 32 indicate the loss in concentration and radiant flux in going from a 
paraboloid concentrator to a Fresnel concentrator with the same aperture 
ratio. Again, the K = 0.01 cylindrical receiver is used. Figure 33 gives a 
description of the Fresnel concentrator. 

If process heating requirements fall within the range achievable by 
Fresnel configurations a foldable Fresnel configuration may be attractive. 

5.5 DESIGN OF A SOLAR CONCENTRATOR TO GIVE A PRESCRIBED 
CONCENTRATION DISTRIBUTION 

As mentioned in Section 3.8 the analytical tools developed can be used 
to design concentrators as well as determine the irradi^tive distribution on 
targets for prescribed concentrator-. To illustrate, a concent; .tor suru 
has been determined which will give a concentration of about 4000 on a section 
of surface of a cylinder of radius R = 0.01. Due to the problems discussed 
previously (in the development of the analysis) it is not always possible to 
develop a design that will give an exact irradiative distribution if an extended 
source is being used. Figure 34 gives the concentrator geometry and Fig. 35 
gives the concentration distribution on the cylinder. The waviness in the dis- 
tr ; bution is due to several factors. The source characteristics and the points 
chosen on the receiver for the numerical calculations are probably the most 
important factors. By manipulating various parameters that must be chosen, 
such as receiver points for the design, the distribution could be made more 
unifoxm. However, the point is that concentrators to give prescribed irradia- 
tive distributions can be designed. 

In Fig. 34 i? e concentrator profile appears to be that of a parabo’-’d. 
Actually the surface is very close to being a paraboloid. Inis is due to the 
fact that the region over which the radiati a is being directed is very small 
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and Fresnel Concentrators 
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relative to the concentrator dimensions. However, the surface in Fig. 34 is 
not a paraboloid as shown in the lower portion of the figure. The quantity 
As is the variation of the surface from a paraboloid, i.e.. 


As 


paraboloid 


s 


5.6 ESTIMATIOt OF CONCENTRATOR SIZES FOR OBTAINING 
MOLTEN ZONES 

The analysis in Section 4 has been used to estimate the power required 
to maintain molten zonet in various size bars of tungsten, silicon and germ- 
anium. With this calculation the solar concentrator that will give the required 
power can be determined* Before performing these calculations a number of 
assumptions must be made and data on the thermal and optical properties of 
the materials obtained. The thermal properties given in Table 3 for the three 
materials can be found in a number of references. However, material surface 
properties such as absorptivity are rather difficult to obtain and the available 
data is rather scattered. Upon examing the data compiled in [64] it was decided 
to assume that a = c = 0.5 for all three materials. 


It aas assumed that the material to be processed was a cylinder and the 
molten zone was two cylinder diameters long. It was also assumed that one 
diameter of length on each end of the molten zone was not insulated. Since 
access to the focal region from the center of the concentrator would perhaps 
be desirable in a design, it was assumed that the inner diameter of the con- 
centrator is one meter. Also, to take into consideration the reflectively of 
the concentrator and any obscuration that might result from process supports, 
etc., it was assumed that >'j = 0.9. 

Except for the region of the molten zone be ig irradiated the process 
material will be radiating to the blackness of space either directly or through 
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Table 3 

MATERIALS PROPERTIES 


i Material -► J 
t Property ■ 

Si 
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k( cal \ 1 

\sec cm °K/ j 
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0.23 | 
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0.09 | 

j 

—i 

' fe) i 
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5 49 1 
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0.5 1 

c 

| 0.5 

j 
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the reflecting concentrator surface. In view of this the radiation geometry 
factors were assumed to be unity. 

With the above assumptions the curves in Figs. 36 and 37 were generated 
using the analysis in Section 4. For the situation chosen the concentrator 
aperture is directly proportional to the bar diameter for a given material. If 
a bar of material 10cm in diameter (about 4 inches) is to be melted we can see 
that about 8kW would be required for silicon, about 18 kW for germanium and 
about 40 kW for tungsten. This corresponds to solar concentrator diameters 
of about 4 meters, 6.5 meters and 29 meters, respectively (for the assumptions 
that have been made). The power requirements and thus the concentrator diam 
eter increases exponentially with melt temperature. This of course is due to 
the radiation losses. 

5.7 SPECIFIC POWER COMPARISONS 

For the space processing of materials on a production basis thermal 
energy for processing must come from the solar source or a nuclear reactor. 

If the solar source is used, the only competitor currently with solar concen- 
trators is solar cells. A meaningful comparison can be made only if complete 
systems are designed. However, the surfaces used to collect the solar radia- 
tion are a major portion of the process system and specific power for these 
surfaces is a good indication of system specific power. The compari* ns 
shown in Table 4 are nwdc for the collecting surfaces. 


In Table 4 specific power comparisons are made for a solar concentrator 
with efficit ncy of 25, 50 and 100%. Based on previous studies, materials 
optical properties, etc., an overall eft'<~iency of 50% is realistic and perhaps 
a little conservative. This 50% efficiency is the actual thermal energy that 
could be used. In the case of solar colls the electrical energy would have to 
be converted to thermal energy and directed onto the melt zone. There would, 
of course, be losses in converting and applying this energy. In view of this the 
solar concentrator is actually more attractive than is indicated in the table. 
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Process Bar Diameter versus Concentrator Aperture 
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Table 4 

SPECIFIC POWER COMPARISONS 


Collecting 

Surface 

Specific Power 
( W att s / Pound ) 

jV 

Solar Concentrator 
Efficiency = 100% 

144.5 

Solar Concentrator' 
Efficiency = 50% 

72.2 

Solar Concentrator 
Efficiency =25% 

36.1 

Solar Cells 

lOkW Panel for Sky lab 

2.63 

** 

Solar Cells 

14 W Unit for ATM Module 

6.8 


2 

Based on 0.87 Ib/ft of projected area in the sun 
direction [ 44] . 

Data obtained from Electronics and Control Laboratory, 
Marshall Space Flight Center. 
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Section 6 

CONCLUSIONS AND RECOMMENDATIONS 

Previous studies that are applicable to the space processing of materials 
with a solar concentrator have been reviewed. Also, an analysis and associated 
computer code has been developed which can be used to determine the heating 
distribution on any surface placed in the focal region of a solar concentrator. 
The computer code can consider optical errors, target placement errors, limb 
darkening effects, etc. This computer code can also be used to design solar 
concentrators to give prescribed heating on a target. A simplified analysis 
is also developed to determine the relation between concentrator size and pro- 
cess size. Based on these analyses and the literature reviewed the following 
observations are apparent. 

• The proper thermal environment is very important 
to process success. 

• Solar concentrators can provide thermal energy for 
space processing* 

• Solar concentrators can be designed to give specific 
thermal environments using the analysis in this report. 

• Solar concentrators can probably give mor^ accurate 
heating control than other heating techniques. 

• The heating distribution on targets in the fo* region 
of solar concentrators is sensitive to target placement. 

• The technology is available for producing large orbiting 
solar concentrators. 

• The technology for pointing solar concentrators accu- 
rately is available. 

• The technology is available for controlling the solar 
concentrator — processing system. 

Although many points that are jnportant to solar concentrator design 
have been covered a detailed design of a solar concentrator system is beyond 
the scope of this study. Also, before a detailed design is arr' ' r ' " l ken firm 
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decisions must be made concerning the material(s) to be processed, size of 
the process, etc. 

The most important conclusions one can draw from this study are that 
solar concentrators are very attractive for space processing, and they cur- 
rently appear more attractive than alternate systems. The first conclusion 
is the result of showing that solar concentrators can provide sufficient 
thermal energy and showing that a great deal of flexibility is available for 
tailoring the process thermal environment. The second conclusion is based 
on the specific power comparison that was made with solar cells. 

From the insight gained from previous studies and the analysis devel- 
oped, the following observations have been made concerning the design of 
solar concentrators for processes that require specific thermal environments. 

• A solar concentrator can be designed to provide the 
thermal environment for a specific process geometry. 

• Multiple -purpose solar concentrators can probably be 
achieved through the use of a changeable secondary 
element (such as the "flux trap 1 ’ mentioned earlier) or 
by using a segmented concentrator whose elements can 
translate along the optical axis. 

If the space processing of materials with a solar concentrator is to be 
a reality, it is recommended that the following efforts be initiated. 


• Terrestrial experimental studies of crystal growth 
and zone refining with a solar concentrator. 

• Develop a detailed analytical model of the process 
heat transfer considering convection, conduction, 
and radiation. 

• An economic and weight comparison of a complete 
solar concentrator-space processing system with 
alternative space processing systems. 

• Initiate a program to provide a space processing 
system and facility for the shuttle program. 
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Appendix A 

A computer code has been developed based on the analyses in this 
report. The program has a great deal of flexibility in that it can determine 
solar concentrator geometries that give a desired heating distribution as well 
as determine the heating on any target placed in the focal region of a concen- 
trator. The highlights of this program are summarized below: 

• Computation of heat flux distribu 4 ion on any receiver 
from specified reflectors 

• Design of continuous or faceted reflectors to give pre- 
scribed heat flux 

• Computation of heat flux distribution on receiver per- 
turbed from desired position 

• Statistical perturbation of ideal reflectors to simulate 
real reflector imperfections, and 

• Limr darkening effects. 


A simplified flow chart, showing the interaction of the main components 
of the program, is presented in Chart A-l. 
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Chart A- 1 - Simplified Flow Chart for SCAP 













